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The gyrokinetic turbulence code GS2 was used to investigate the effects of plasma f3 on linear, coUisionless 
ion temperature gradient (ITG) modes and trapped electron modes (TEM) in National Compact Stellarator 
Experiment (NCSX) geometry. Plasma /3 affects stability in two ways: through the equilibrium and through 
magnetic fluctuations. The first was studied here by comparing ITG and TEM stability in two NCSX 
equilibria of differing /3 values, revealing that the high /3 equilibrium was marginally more stable than the low 
/3 equilibrium in the adiabatic-electron ITG mode case. However, the high /3 case had a lower kinetic-electron 
ITG mode critical gradient. Electrostatic and electromagnetic ITG and TEM mode growth rate dependencies 
on temperature gradient and density gradient were qualitatively similar. The second /? effect is demonstrated 
via electromagnetic ITG growth rates' dependency on GS2's /3 input parameter. A linear benchmark with 
gyrokinetic codes GENE and GKV-X is also presented. 



I. INTRODUCTION 



Magnetic fusion energy requires the containment of 
very hot plasmas for a long enough time to allow fu- 
sion reactions to occur. Turbulent transport (most likely 
the result of drift-wave instabilities) breaks this confine- 
ment and can cause a significant amount of heat loss 
in tokamaks and spherical toriii In contrast, neoclassi- 
cal transport can often account for the poor confinement 
in traditional stellar atorSi^ However, modern stellarator 
designs, such as Wendelstein 7- AS (W7-AS),^ Wendel- 
stein 7-X {Wl-X),"^'^ the National Compact Stellarator 
Experiment (NCSX) ,6 the Large Helical Device (LHD),i 
and the Helically Symmetric Experiment (HSX)^"— have 
shown or are designed to have improved neoclassical con- 
finement and stability properties. Thus, plasma tur- 
bulence and transport levels may be experimentally- 
relevant now and could affect performance of these stel- 
larators. 

Gyrokinetic studies of drift-wave-driven turbulence in 
stellarator geometry are relatively recent and compre- 
hensive scans are scarce. Most of these studies were done 
using upgraded versions of well-established axisymmetric 
codes: the linear eigenvalue FULL code^itiiS, the nonlin- 
ear initial- value or eig envalue GENE code^ii^ and the 
nonlinear initial-value code used in this paper, GS2. The 
nonlinear initial-value GKV-X code^^iii which uses the 
adiabatic electron approximation, was specifically writ- 
ten to simulate turbulence in stellarator geometry. All 
four codes use the flux-tube limit in their geometry, al- 
though GENE has been upgraded to allow for full flux- 
surface simulations. The microinstability code GS2^^ was 
extended from its original axisymmetric-geometry ver- 
sion to treat the more general case of non-axisymmetric 
stellarator geometry as described by Ref. il3and briefly 
mentioned in Section [III The work in Ref. Il9| includes 
a linear benchmark of GS2 stellarator simulations with 



FULL. Section [ITT] briefly displays the results of a linear 
benchmark of GS2 with GENE and GKV-X, rounding 
out the major stellarator codes. 

Finally, in Section |IV] the upgraded GS2 is used for 
comprehensive parameter scans and instability studies in 
the National Compact Stellarator Experiment (NCSX) 
design. NCSX, with its quasi-axisymmetric magnetic 
conflguration, is a bridge in configuration space between 
tokamaks and the rest of the stellarator world. Therefore, 
it is an excellent configuration to begin detailed gyroki- 
netic stellarator studies with GS2, which has been used 
successfully on axisymmetric geometry for many years. 
The two NCSX configurations used in Section IIVI were 
created as part of a series of flexibility studies^" that were 
performed using a magnetic coil set similar to the flnal 
design of the NCSX machine. The equilibrium optimiza- 
tion code STELLOPT was used to find currents in these 
coils needed to meet desired configuration properties. As 
a consequence, sets of configurations exist in which only 
one parameter, such as magnetic shear and plasma /3, 
varies significantly. Studies in Section IIVI survey linear 
stability in two configurations that differ only by plasma 
equilibrium /3. These equilibria were compared via the 
growth rates of the electrostatic adiabatic ITG mode, 
electrostatic coUisionless kinetic ITG-TEM mode, and 
electromagnetic coUisionless kinetic ITG-TEM mode. 



II. GS2 COORDINATE SYSTEM 

First, GS2 geometry input must be built by a series 
of programs. VMEC^^'^^S creates the 3D MHD equilib- 
ria used as a basis for all gyrokinetic stellarator codes. 
Qjgrp23 extracts from the full 3D equilibrium the geo- 
metrical data needed to represent a fiux tube based on 
field-line-following coordinates. This coordinate system 
includes the radial coordinate, p = y/s (s ss (r/a)'^ is the 
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normalized toroidal flux), the distance along a field line, 
9, and the angle that selects a flux tube, a = C,—q{0 — 9Q) 
(where C and are the Boozer toroidal and Boozer 
poloidal coordinates and 9^ is the ballooning parameter) . 
To obtain the final GS2 geometry input file, FIGG^ uses 
the GIST output file to calculate the pitch angle parame- 
ter grid, involved in the velocity integration of the distri- 
bution function. (The pitch angle parameter, A = /^/-E, 
is related to v\\ through v\\/v — \/l — XB.) 

The GS2 documentation^^ defines geometrical quanti- 
ties in terms of a parameter n / dp, where p is the radial 
coordinate and at is the normalized poloidal fiux. Geo- 
metrical quantities in this paper follow GS2 notation and 
include d'^N /dp. For more information, see Refs. ^ and 



III. BENCHMARKS OF GS2, GENE, GKV-X 

GS2, GENE, and GKV-X results were compared for an 
NCSX VMEC equilibrium based on the standard S3 con- 
figuration of NCSX design. This configuration is quasi- 
axisymmetric with three field periods, an aspect ratio 
of 3.5, and a major radius of 1.4 m. The following 
benchmark used geometry with the surface at s = 0.5 
(r/a ~ 0.7), the a = field line, and the ballooning pa- 
rameter ^0 = 0. The average /3 is = 4% and, at this 
surface, the safety factor is q = 1.978. 

Figure [T] shows the variation of the magnitude of the 
magnetic field along the chosen magnetic field line, with a 
resolution of 209 9 grid points per poloidal period. There 
were approximately 30 points in the pitch angle param- 
eter grid. The 9 range extends from — 37r to 37r. 

Figures [2] and [3] are the variations of {k±/n)'^, where 
n is the toroidal mode number, and the curvature drift 
along the same chosen field line. By convention, positive 
curvature drifts are bad or destabilizing, while negative 
curvature drifts are good or stabilizing. Significant un- 
stable modes occur where k± is small, which is near 9 = 
for this equilibrium, since instabilities are generally sup- 
pressed at large k± by FLR averaging. Because Figure 
[5] indicates that the curvature is bad in this region near 
9 = 0, where k± is the smallest, it is expected that un- 
stable modes will appear here. 

The benchmark case was an electrostatic, collision- 
less ITG mode with adiabatic electrons. The temper- 
atures were such that Ti = Tp, the temperature gradient 
was ajqjLT = 3 (ajqlLx = —aN{^IX)dX/dp), and the 
density gradient was / = 0, where the normaliz- 
ing length was chosen to be an averaged minor radius, 
aN~ 0.323m. See Table HI 

Figure |3] shows the growth rate and real frequency 
spectra for this mode. The maximum discrepancy in 
growth rate between GS2 and GKV-X is 8%, with GENE 
always in between. The frequencies agree to within 5%. 
This agreement is excellent. GS2 and GKV-X's electro- 
static potential, 4>^ for a particular kyPi = 0.9 is shown 
in Figure [S] These electrostatic eigenfunctions also agree 
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FIG. 1. B vs. e grid for NCSX QAS3, with s = 0.5, a = Q, 
and 9o — 0. (color online) 




e (rad) 

FIG. 2. Variation of the curvature drift term (lJcv = 
{d-^N/dp){\i±/n) ■ b X [b ■ Vb]) (for n = 1) along 9 for NCSX 
QAS3, with s = 0.5, a — 0, and 9o = 0. (color online) 



well. 

Now, all four gyrokinetic stellarator codes have been 
benchmarked linearly (here GS2, GENE, and GKV-X; 
Ref. 19 describes GS2 and FULL's benchmark). In ad- 
dition to those in NCSX geometry, comparisons of linear 
GENE and GS2 resuhs for W7-X also agreed well.^^ With 
these successful benchmarks, further studies can be per- 
formed with more confidence. The next section begins 
this venture with GS2. 



IV. NCSX /3 STUDIES 

High plasma /3 is important for fusion because the fu- 
sion power at fixed magnetic field is approximately pro- 
portional to /3^. To start studying the effect of plasma 
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FIG. 3. Variation of (^)2(6i)(d*jv/dp)2 for NCSX QAS3, 
with s — 0.5, a = 0, and Oq — 0. (color online) 
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Ro 


~ 4ajv ~ 1.4m 


aN 


^ 0.323m 


Ba = {B) 


^ 1.6T 


rrlref 


2mp 


Vt 


y^Ti/rriref 


GS2 cj units Vt/a,N 


^ 6.782 X lO'^sec"^ 



TABLE I. The set of local parameters used the microinstabil- 
ity simulation based on the NCSX QAS3 equilibrium. 



beta on gyrokinetic turbulence in stellarators, linear ITG 
and TEM stability was compared for two configurations, 
one with equilibrium /3 = 0% and one with = 4%. 
Figure 6 of Ref. shows the poloidal cross-sections for 
three toroidal locations, along with t profiles for vari- 
ous plasma currents (Ip). The plasma shape varies very 
little with Ip. This section uses set of beta scans with 
Ip = 17AkA, because their t profiles varied the least, al- 
lowing for isolation of the effects of /3. 



A. Discussion of GS2 (3 Parameter 

The physical beta enters into GS2's equations (the 
gyrokinetic and Maxwell's equations, see Ref. [26l ) in 
two main ways, through its indirect effect on the MHD 
equilibrium (such as the Shafranov shift and the cur- 
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FIG. 4. Variation of 7 and ujr with kypi for NCSX QAS3, 
comparing three codes. 




FIG. 5. Comparison of vs. 9 (with kyPi = 0.9) for GS2 
iRe{(j)): blue line, Im{(j)): black line) and GKV-X {Re{(j>): 
green circles, Im{(p): red crosses), (color online) 



vature drift) and through its direct effect in the gy- 
rokinetic equations, controlled through the parameter 
Pinput- This GS2 beta parameter is defined as Pinput = 
2fionrefTref / B^^j:, the ratio of the reference pressure to 
the reference magnetic energy density. Pmput is used 
in the scaling of SB^^ = V_l x A_l and SA^^, through, 
for example, the weighting of the contribution of each 
species to th e total parallel current by a factor Ws = 
'^PinputZsrisyjTs/ms. While Pinput must be set to match 
Pequii in the geometry files for consistent results, setting 
Pinput = is a convenient way to turn off magnetic fluc- 
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FIG. 6. NCSX beta flexibility studies comparing \B\ vs. for 
both P = 0% and /3 = 4%, at ^ = r/a = 0.5, a = 0, and 
^0 = 0. (color online) 



FIG. 7. NCSX beta flexibility studies comparing the curva- 
ture drift terms {oJcv,norm = {2a% / BN)id'^ n /dp){k±/n) -h x 
[h ■ Vb]) along 6, for 13 = 0% and 4%, at ^ = r/a = 0.5, 
a = 0, and do = 0. (color online) 



tuations and focus only on electrostatic fluctuations, as 
done in Sections IIVCWdI 



B. Geometry and Plasma Parameters 

For these (3 studies, the geometry used had the surface 
with normalized toroidal flux s ~ ((r/a))^ — 0.26, field 
line a = 0, and ballooning parameter 6o = 0. The mag- 
nitude of magnetic field, curvature and V-B drift com- 
ponents, and \k±\'^ along the field line are plotted for 
both (3 = 0% and /3 = 4% in Figures EHSl More param- 
eters for both equilibria are in Table HH All growth rate 
and frequency values are normalized such that (7,w) = 
{^'fphy steal , ^physical ){a/vthi)- These runs are coUisionless 
(collision frequency v = 0). For the following studies, 
several plasma parameters were varied around the base 
case parameters shown in Table IIIII 

For each equilibrium (/? — 0% and (3 ~ 4%), conver- 
gence studies were run with increasing resolution in 9 and 
velocity-space for single ion species, ITG-driven adiabatic 
electron modes. Here, GS2 studies use grids with 30 pitch 
angle parameter points for both the (3 — 0% and 4% 
equilibria, approximately 750 9 points for /3 = 4%, and 
approximately 630 9 points for /3 = 0%. These grids' 
results were well- converged (to within a few percent of 
the results from higher- resolution grids). There are 32 
energy grid points. 



Electrostatic Adiabatic ITG mode 



Using the base parameters in Section |IVB[ linear ITG 
stability as a function of temperature gradient, un/Lt, 
was compared for both equilibria, over a wavenumber 
range of kyPi G [0.6,1.4] (note in Figure IH] that this 



150 



100 



50 








-10 



-5 





e (rad) 



10 



FIG. 8. NCSX beta flexibility studies comparing (^^'J vs. 

e for /3 = 0% and 4%, at = r/a = 0.5, a = 0, and Oq = 0. 
(color online) 
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2.175 


2.011 
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0.356 


0.278 


W) 


0.0% 


4% 


R 


~ 4.7ajv ~ 1.5m 


~ 4.7a]v ~ 1.5m 


ajv 


X 0.322m 


^ 0.322m 


Ba = (B) 


i.58r 


i.55r 



TABLE II. Geometry values for the NCSX /3 = 0% and 4% 
equilibria. 
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G [0.6, 1.4] 


T, = T, 


IkeV 


V 





rriref 


2mp 


Vt 


^yTi/rriref 


GS2 cj units vt/aN 


^ 6.214 X lO^sec"^ 



TABLE III. The base set of local parameters used in the 
NCSX P studies. 



range is sufScient to capture the peak of the growth 
rate spectrum, and therefore the fastest growing mode). 
The peak growth rates for both equihbria occur between 
kyPi « 0.6 and « 1.2 and are shown in Figure [Til indicat- 
ing that the critical gradient of the /? = 0% equihbrium 
is dN / LT,crit ~ 1.13 and that of the f5 — 4% equihbrium 
is ttN / Lr.crit ~ 1.16. These values are not significantly 
different. The fact that beyond marginal stability, the 
growth rates of /3 = 0% are larger than those oi fi — 4% 
is an indication that /? is stabilizing to the ITG mode, as 
found in the tokamak studies in Ref. A representa- 
tive eigenfunction is shown in Figure [TUJ note the typical 
ballooning-about-zero behavior of an ITG mode. 

Looking at the effect of the density gradient on the 
critical temperature gradient in Figures [T^llTSl with 
0.N I Ln = 1, CLN / Lx^crit lowcrs by w 0.1 with respect 
to the aN / Ln = value in each case, appearing to be 
somewhat destabilizing. aN / Ln > 2, however, appears 
to be strongly stabilizing, consistent with a transition to 
the slab limit of the ITG mode where a density gradient 
is stabilizing.^^ 

As expected, when one compares the growth rates as a 
hmction of aj^ / for various values of a-^lLj- (Figs. 
[rilfTS]) . the growth rates increase monotonically with 
um/Lt- Also, the growth rates for (3 = 0% are higher 
than those for /3 = 4%, another sign that higher (3 could 
be somewhat stabilizing. 



D. Electrostatic Kinetic ITG-TEM 
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FIG. 9. ITG adiabatic electron growth rates vs. kypi for 
NCSX /3eq«« = 0%, fliv/I/n = 0, and ajv/LT = 1.6. (color 
onhne) 
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FIG. 10. A representative electrostatic potential eigenfunc- 
tion for these ITG modes. ajv/I/T = 1-6, aN/Ln = 0, and 

kypi — 1. 



Adding kinetically-treated electrons allows one to 
study the trapped electron mode (TEM) and hybrid ITG- 
TEM (driven by both / Lt and a]sr/Ln) modes. Figure 
[TBI shows the kyPi spectrum for ajq/Ln = I^un/Lt = 
0.5 and Figure [17| the kyPi spectrum for aN/Ln — 
2,aN/LT = 0. The peak of the growth rate spectrum 
shifts as aN/Ln and aN/L^ change. When the gradi- 
ents are large enough (Figure [T7| , two distinct regimes 
are seen, with a mode switch evident in the frequencies. 
These studies focus on the peaks lower than kypi = 1.8. 

Figures [T51[T^ show growth rates vs. gn/Lt (where 
aN / Lt — qn/Ltz — qn/Ltc) for several values of 
o-N/Ln, for both the equilibrium with j3 = 0% and that 
with /? = 4%, for the kyp, G [0.4,1.8] with the highest 
growth rate. Both have the same general trend: for all 
values of qn / Ln, past a critical temperature gradient, the 



growth rates increase almost linearly with un/Lt, indi- 
cating that this mode is driven by the temperature gra- 
dient. When aAr/L„ = 0, in both cases, there appears to 
be a critical temperature gradient, which is lower than in 
the adiabatic electron case. Here, aN / LT,crit.0=o ~ 0.75 
and ttN / LT^crit,/3=4: ~ 0.25. Also as in the adiabatic case, 
increasing aN/Ln first further destabilizes the mode-the 
large linear growth begins for a lower temperature gra- 
dient than for aN/Ln = 0-but then it is stabilizing for 
higher density gradients (this is more easily seen in Fig- 
ures [20H1I]) ■ Though, as density gradient increases, the 
value of the flat part of the growth rate for low qn/Lt 
increases: this is a density-gradient-driven regime. Com- 
paring the two /3 equilibria, the /3 = 4% growth rates 
for the fiat part of the plot are lower than the /3 = 0% 
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FIG. 11. ITG adiabatic electron growth rates vs. temperature 
gradient for NCSX Pequu ~ 0% (blue circles) and 4% (red 
crosses), aN/L„ = 0. Fits (dashed lines) obtained through 
piecewise linear interpolation on the lowest half of the growth 
rate curve, (color online) 



FIG. 13. Electrostatic ITG adiabatic electron growth rates 
(at kyPi of maximum 7) vs. temperature gradient for NCSX 
Pequii = 4% for various density gradients: aN / Ln = (blue 
circles), at^jLn = 1 (green dashed line), aj^jLn = 2 (red 
crosses), (color online) 





FIG. 12. Electrostatic ITG adiabatic electron growth rates 
(at kypi of maximum 7) vs. temperature gradient for NCSX 
Pequii = 0% for various density gradients: un / Ln = (blue 
circles), aN / Ln ~ 1 (green dashed line), a^ILn ~ 2 (red 
crosses), (color online) 



FIG. 14. Electrostatic ITG adiabatic electron growth rates (at 
kyPi of maximum 7) vs. density gradient for NCSX P^quu ~ 
0% for various temperature gradients: un/Lt ~ 1.6 (blue 
circles), un/Lt ~ 2.0 (green dashed line), on/Lt ~ 2.4 (red 
crosses), gn/Lt ~ 2.8 (black triangles), (color online) 



case, for a^f /L^ — 0, 1, and higher for apf/Ln — 2. It ap- 
pears that the "critical gradients," for the strong linear 
growth at higher un/Lt, are lower, as well as the critical 
gradient for aN /Ln — 0. 

These results seem to differ some from Ref. [2^, which 
appears to find a larger region of stability for sufHciently 
small aN / Ln and aN/LT, for a particular set of tokamak 
parameters. Ref. [2^ however, included finite collisions, 
while these studies are collisionless. Including finite colli- 
sions in these simulations may increase the stability win- 
dow for TEM at weak gradients, as has been found in 
tokamaks^o and STs^ 



Figures HOIHT] show growth rates vs. aN /Ln for several 
values of aN/LT (again for both equilibria and for the 
kyPi e [0.8, 1.4] with the highest growth rate). Here, for 
gn/Lt > 1, one can more easily see the increased desta- 
bilization of the mode as aN/Ln increases, until about 
o-N/Ln — 1, when the growth rate decreases. For val- 
ues of the temperature gradient lower than the adiabatic 
electron critical temperature gradient of oat/Lt ~ 1-3, 
the mode is density-gradient driven: the growth rate in- 
creases slowly with aAr/L„. The growth rates are again 
higher for the /3 = case. 
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FIG. 15. Electrostatic ITG adiabatic electron growth rates (at 
kypi of maximum 7) vs. density gradient for NCSX Pequu = 
4% for various temperature gradients: un/Lt = 1.6 (blue 
circles), aN/Lr = 2.0 (green dashed line), aN/Lr = 2.4 (red 
crosses), aN /Lt = 2.8 (black triangles), (color online) 



FIG. 17. Electrostatic ITG-TEM kinetic electron growth 
rates vs. kypi for NCSX Pequu ~ 0%, aN / Ln = 2, an / Lt = 0. 
(color online) 
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rates vs. temperature gradient for NCSX Pequu ~ 0% for var- 
FIG. 16. Electrostatic ITG-TEM kinetic electron growth ious density gradients: ajv/L„ = (blue circles), a]v/L„ = 1 
rates vs. kyPi for NCSX jS^q^a = 0%, aN/L„ = 1, qn/Lt = (green dashed line), ajv/L„ = 2 (red crosses), (color online) 
0.5. (color online) 



E. Electromagnetic simulations 

As a preliminary investigation of electromagnetic ef- 
fects, the GS2 beta parameter, Pinput, was scaled using a 
fixed equilibrium, with temperature gradient a^q / Lt — 5. 
For all following discussions, the notation used is Pgs2 = 
'^Pinput, to convert from /? for a single species to per- 
cent /? for two species (an electron and an ion species of 
equal T and n). In order to demonstrate the effect that 
instability-driven current fluctuations have on the growth 
rate, equilibrium /3 was held fixed and Pgs2 scanned. Fig- 
ure [21] compares two /3gS2 scans based on configurations 
with equilibrium /3 = 0% and 4%. The frequencies and 



growth rates match closely when PgS2 — 0%. But, the 
apparent mode switch occurs earlier for the /? — 0% equi- 
librium (around /3gS2 = 1.5%) than the /? = 4% equilib- 
rium (/3g52 — 2.0%). In addition, the values in growth 
rate and frequency differ by about 20% when /3gs2 = 4%, 
indicating that matching this GS2 parameter with the 
equilibrium value does matter. The general trend, similar 
to tokamak results, is that (5gs2 is stabilizing to the ITG 
mode at moderate values, but the fastest growing mode 
switches character to a high frequency mode (perhaps a 
kinetic ballooning mode) at higher PgS2- Equilibrium /3 
is stabilizing for this higher frequency instability (this is 
the stabilizing mechanism that can give rise to the second 
stability regime for MHD ballooning modest). 
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FIG. 19. Electrostatic ITG-TEM kinetic electron growth 
rates vs. temperature gradient for NCSX P^q^u = 4% for var- 
ious density gradients: an/Ln = (blue circles), aN / Ln — 1 
(green dashed line), aN/Ln = 2 (red crosses), (color online) 
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FIG. 20. Electrostatic ITG-TEM kinetic electron growth 
rates vs. density gradient for NCSX PequU = 0% for various 
temperature gradients: un/Lt ~ (blue circles), un/Lt ~ 2 
(green dashed line), um/Lt = 3 (red crosses), (color online) 



Figures [23ll24l show the magnetic potential, A^, for 
f!^equii = 4% and /3gs2 = 1% and 4%. This demonstrates 
the effect PgS2 has on the perturbed magnetic fields-the 
magnitude of A|| is much larger in the /3gS2 = 4% case 
than in the /3gS2 = 1% case. Figure [25] is the magnetic 
potential, A\\, for f3equU = 0% and /3gS2 = 4%. Note 
that it is almost identical to Figure [Ml l^equU — 4% and 
PgS2 = 4%, demonstrating that only the /3gS2 parameter 
affects the fluctuating magnetic fields, not PequU- 

In electromagnetic GS2 runs, one always includes 
5Bi^ = VA|| X z when calculating the fields, but one can 
choose to include 5B\\ = V_l x A_l— or set it to zero. One 
might want to ignore this term to save computational 
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FIG. 21. Electrostatic ITG-TEM kinetic electron growth 
rates vs. density gradient for NCSX jigquu ~ 4% for various 
temperature gradients: aN / Lt = (blue circles), aN /Lt ~ 2 
(green dashed line), aN /Lt = 3 (red crosses), (color online) 
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FIG. 22. Growth rates and real frequencies, in normalized 
units (7GS2,t^GS2) = {■y,Uj)/(vT /a) as a function of ^gs2, for 
equilibrium /? = 0% (uj: red solid line, 7: blue dashed line) 
and P — 4% {uj: blue circles, 7: green crosses). kyPi — 1.0. 
(color online) 



time. Figures [26l[27l demonstrate the importance of in- 
cluding (5i3|| for high /3gs2 values. For /?gs2 ^ 1.5%, the 
growth rates and frequencies for SB^^ — and SB^^ ^ 
are approximately equal, because /3gS2 scales the SB\\ 
field, so that when /3gS2 is low, the contribution from SB\ 
is small. However, as PgS2 increases past Pgs2 = 2%, 
the contribution from (5_B|| increases: including (5_B|| has 
a destabilizing effect at higher Pgs2- 
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FIG. 23. A|| for = 4%, /3gs2 = 1%, kyp, = 1.0. Blue: 

7m(j4||), Green: Re{A\\). (color online) 



FIG. 25. 4|| for = 0%, /3gs2 = 4%, 

Im(j4||), Green: iie(^||). (color online) 
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FIG. 26. NCSX Pequii = 0%: growth rates and frequencies 
((7GS2,tJGS2) = {'y,oj)/{vT /a)) vs. Pgs2 for / (cj: blue 
circles, 7: green crosses) and 5B\\ = (u: red solid line, 7: 
blue dashed line), (color online) 



CONCLUSION 



Understanding the effects of stellarator geometry on 
gyrokinetic turbulence, along with how much of an effect 
turbulence has on confinement in current experiments, is 
of paramount importance for designing future magnetic 
condiment fusion devices. Drift- wave instabilities are be- 
lieved to cause turbulence, and can be modeled using sev- 
eral gyrokinetic codes, including GS2. The nonlinear gy- 
rokinetic turbulence code GS2's non-axisymmetric geom- 
etry capabilities were linearly benchmarkcd for an NCSX 
equilibrium with GENE and GKV-X. The growth rates 
and real frequencies of an adiabatic ITG mode agreed to 
within 8% for all three codes. Coupled with a previous 
GS2 benchmark with FULLfiS, all four gyrokinetic stel- 
larator codes have been benchmarked successfully against 
each other. 



Extensive studies of instabilities in two NCSX equi- 
libria were conducted. Comparing NCSX equilibria of 
differing j3 values revealed that the /3 = 4% equilibrium 
was marginally more stable than the /3 = 0% equilib- 
rium in the adiabatic-electron ITG mode case, but less 
stable in kinetic electron ITG-TEM mode case. How- 
ever, their electrostatic adiabatic ITG mode and electro- 
static kinetic ITG-TEM mode growth rate dependencies 
on a^v/iT, cln/Lu, and kyPi were similar. 

There are two effects of finite plasma /3 on microin- 
stabilities. The first is created by the changes in mag- 
netic geometry, and this affects even electrostatic modes. 
The second effect is due to fluctuating currents (and 
magnetic fields), and this can be varied as an indepen- 
dent parameter in calculations (although not in the real 
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tently with the equilibrium /3 in order to have physi- 
cal results. For a fixed magnetic equilibrium, the ef- 
fect of /3gs2 on magnetic fluctuations is at first stabi- 
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pabilities of all of the involved gyrokinetic codes. 
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